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The magnetic moment and vanishing charge of a Dirac neutrino are physically observable 



I quantities and must not depend on the choice of gauge in a consistent quantum field theory. We 

pL| ■ verify this statement explicitly at the one loop level in both and unitary gauges of the min- 

I imally extended standard model. We accomplish this by manipulating directly the integrands 

of loop integrals and employing simple algebraic identities and integral relations. Our result 
' r-|'' generally applies for any masses of the relevant particles and unitary neutrino mixing. 
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1 Introduction 



The fundamental properties of a particle like its charge and electromagnetic dipole moments 
are physical quantities that in principle are experimentally measurable, for reviews, see Refs. 
|[I1I3- These quantities can be unambiguously calculated in a quantum-mechanically consistent 
theory like the standard model (SM), and confronted with the measurements to decide whether 
the theory is correct or not. That said, the practical calculation and demonstration of its result 
being independent of computational methods are not always trivial. We have witnessed a similar 
circumstance recently, concerning the one-loop contribution of the charged weak gauge bosons 
to the two-photon decay rate of the Higgs boson. A new computation in unitary gauge [l3l|4l 
claimed an answer that is different from the well-spread result obtained long ago JH |6l |7l [51 in 
a special gauge, i.e., the 't Hooft-Feynman gauge = 1) among the class of renormalizable 
gauges. Subsequent studies by various methods, including computing in both and unitary 
gauges, see for instance Refs. [|9l [TOl, confirmed the old result, and taught us a great deal on 
computational subtleties in a theory that is nontrivial in the high energy regime. 

In this work we examine a similar problem in the neutrino sector, i.e., the charge and mag- 
netic moment of a Dirac neutrino in SM that is minimally extended by the introduction of 
right-handed neutrinos. We show explicitly at the one-loop level in both R^ and unitary gauges 
that the neutrino charge vanishes and its magnetic moment is a gauge independent quantity. 
The issue has been partially studied in the literature. Early works [fTTl [T2l [T3l [T4l [T5l assumed a 
massless neutrino or expanded the quantities to the leading order in the small masses of neutri- 
nos and charged leptons, ignored the lepton mixing, or computed in a special gauge. A further 
step was taken some years ago [|T6l[T7llT8l . It was found [fTTl . for instance, that up to the second 
order in the expansion of small neutrino masses the charge vanishes and the magnetic moment 
is (^-independent, and that the charge vanishes exactly in the 't Hooft-Feynman gauge. Here 
we cope directly with the integrands of loop integrals, and demonstrate manifestly that both 
quantities are gauge independent for any masses of the relevant particles and for any unitary 
lepton mixing. 

In the next section we set up our notations and suggest how to calculate in a nice way to 
isolate the terms that potentially contribute to the charge and magnetic moment. We describe 
in some detail in sec [3] our calculation in R^ gauge. Our one-loop exact result for the magnetic 
form factor at vanishing momentum transfer is shown in Eq.(l59l). This is followed by a short 
discussion in sec|4]on the calculation in unitary gauge. We summarize briefly in the last section. 



2 Computational strategy 



The charge and magnetic moment of a Dirac particle can be defined by the amplitude of a 
process in which it radiates a photon. 



u{p^)i£^^{q)u{p^ 



-ie)u{p-) y^Fi{q^)-^io^yq^F2{q^) + --- u{p^ 



(1) 
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Here p± — p±q/2 are the momenta of the incoming and outgoing particle of mass m, and q 
is the photon's outgoing momentum. The above decomposition in terms of the standard form 
factors is based on Lorentz covariance and electromagnetic gauge invariance, and assumes that 
the Dirac particle in both initial and final states is physical: 

i>±u{p±) = mu{p±), p\ = n?. (2) 

The dots in Eq. ([T]) stand for two more form factors that are irrelevant here; one corresponds 
to the electric dipole moment that cannot occur at one loop in the minimally extended SM (as 
can also be seen from sec|3]), and the other is the so-called anapole whose Lorentz structure is 
quadratic in q. The form factors at an arbitrary q^ are generally not measurable quantities, since 
the above (unphysical) amplitude appears as part of the complete contribution to a physical 
process. Nevertheless, Fi(0) and ^2(0) are physical quantities because they correspond to the 
charge and anomalous magnetic moment of the particle. Our convention is such that the electron 
has the charge eFi(O) = e < and the magnetic moment vector, /l = (e/m)[Fi(0) -|-F2(0)]5 
with S being its spin vector, that appears, e.g., in the interaction potential of the dipole with an 
external magnetic field B,V = — /i ■ B. 

The charge, Fi{0), is relatively easy to isolate. Setting ^ = removes all other Lorentz 
structures, and allows us to employ the equations of motion (EoMs) in the limit ^ — )■ for 
both initial and final particles, ^u{p) = mu{p), to reduce the amplitude completely to the 7^ 
form. There are several ways to work out the anomalous magnetic moment F2(0). One could 
isolate by brute force terms contributing to the form factor F2{q^) and take its value at q-^ = 0. 
Most studies in the literature follow this approach. In the second approach, one employs a 
projection operator, and expresses F2(0) as a combination of Dirac traces [[T9l l20l |2T| . Here we 
take a third approach, which might be the best to observe the cancellation of gauge dependence 
among various Feynman graphs. As we will show in the next section, the cancellation happens 
at the level of loop integrands. In this approach, we take the derivative of the amplitude with 
respect to the photon momentum, idy£^^{q), anti symmetrize it in the Lorentz indices /i and v, 
and then evaluate iidX q = 0. Since all form factors are smooth at q^ = 0, only the magnetic 
moment term survives the procedure and yields —e/{2m)o^vF2{0). [We remind once again that 
the electric dipole term vanishes at one loop but would appear at higher orders.] Comparison of 
the two gives the answer for ^2(0) . 

An important point in implementing the above procedure should be noted. We mentioned 
that the decomposition in eq ([T]) is possible only upon using EoMs When computing ^2(0), 
we are essentially expanding £^^{q) in small q and isolating its linear terms. A term that is 
manifestly linear in q cannot avoid our eyes, for which we are free to apply the limiting EoMs, 
]/)u{p) — mu{p), because the difference to the exact ones does not affect ^2(0). With terms of 
apparently zeroth order in q we should be careful. For these terms, when necessary, we must 
apply the exact equations ^ since the difference now is exactly what we are interested in and 
may enter F2(0). Ignoring this will result in an incorrect, gauge-dependent answer. Another 
point is more technical. Although antisymmetrization in Lorentz indices is not mandatory since 
it will come out automatically upon finishing the calculation, one can simplify the algebra by 
doing antisymmtrization at an early stage. 
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Fig. 1 Feynman diagrams contributing at one loop to the vertex function ir^{q). Wavy 
(dashed, dotted, solid) lines stand for the gauge boson (scalar, ghost, fermion) fields. 
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Fig. 2 Feynman diagrams contributing at one loop to the 7Z mixing energy zn^v(^)- 
There are two classes of Feynman graphs in SM that contribute at one loop to the amplitude 
i^^{q), through the proper vertex ir^{q) in Fig. 1 and the photon-Z boson mixing energy in 
Fig. 2, m^vi^l)^ attached to the tree level neutrino-Z vertex (see the last graph in Fig. 2). While 
the former contributes to both Fi (0) and ^2(0), the latter contributes only to Fi (0) through 



«n^v(o) 



I lg2 

m| 2cw 



fPi. 



(3) 



Here we use the standard notations of SM: m^ z are the masses of the and Z bosons, g2 is the 
gauge coupling of SU (2)l, cw = cos 9w and sw = sin 9w with 9w being the weak mixing angle, 
and Plr = (1 =F Ts) /2. We display here the contributions from individual graphs. Working in 
J-dimensions, we write 



^n^v(o) 



1 



egl\Va^\^ 
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(4) 



(5) 
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where, denoting k± = k±q/2, from Fig. 1, 

(la) = +y^PL{li + i> + ma)ypPLTafi^{-k-M.~q)P'"'{k-)P^''{k+)P-\ (6) 

,,,, , _2{miPL-maPR){Jj^ + ^ + ma){miPR-maPL){k^+k+)^ 
[lb) = -2 ^ Y\7T2 i ' 

, _ , rcjPL{]^ + l/> + 'na){m,PR-maPL)P^^{k+) 



, imiPL-maPR)i]^+^ + ma)rpPLP^'P{k-) 



. ^ yaPL{]i + i>-+ma)y^{l^ + U + ma)ypPLPP^{k) 

[{k + p^f-ml][{k + p^f-ml] ' ^""^ 

_ , ^2 i^jPL - mgPR) {1^ + ^-+ ma)y^ {^+^+ + ma){miPR - rnqPt) 
" ^'"^ [{k + p^r-ml][{k + p.r-ml]Q2 ' ^ ^ 

and from Fig. 2, 

(2a) = -cl/[2gp(jg^y-gp^g(yy-gpvg(y^]P'^'^, (12) 

(2Z>) = -(4-52,)^^ve2'' (13) 

(2c) = +cirpa^,{-k,k,0)rii^y{-k,k,0)P''l^pP'', (14) 

{2d) = +2{cl - sl)k^kyQ2\ (15) 

(2e) = +2slmlP^yQ2\ (16) 

(2/) = -2c^wk^kyQ2\ (17) 

Note that the fermion loop in Fig. 2(g) is transverse and drops out at ^ = 0. We have defined 
the shortcuts for the propagators and triple-gauge vertex: 

^ali^iiPuP2,P3) = {P2-P3)agp^i + {P3-Pl)l58iia + {Pl -P2)iigap: (18) 
Pfivip) = gfiv[p -my^V -SwPnPv[p -t,wm^]' [p -niy^Y , (19) 
P={k + pf-ml, Qi=k^-ml/, Qi = k^ - ^wml/ , (20) 

with P^p = P„p (k) and 5w = l—^w- i^a (mi) is the mass of the charged lepton £a (neutrino v,), 
and Vai is the lepton mixing matrix appearing in the charged current interaction. A summation 
over all £a is always implied. The identical initial and final neutrino satisfies EoMs Q where 
now m = m,. The above loop integrands will be manipulated in the next two sections. 



3 Evaluation in gauge 
3.1 Charge 

Let us start with the charge. Setting ^ = simplifies significantly the expressions of (Ix). 
Using {]^ + ^)r^i{H^) = -{k + p)^r^i+Ak + p)^{li^ + ^), d^Q2^ = -2k^Qf, and d^p-^ = 
—2{k + p)^P^^, {\b) and (1/) sum to a total derivative: 

[(lZ7) + (l/)]o = -my,^d^{[{^+]/>){m}PR + mlPi:)-mml]{PQ2)-']. (21) 
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where the subscript denotes evaluation at q = 0. Considering the relation 

rap^{-kXO)P"PP^'' = {kPp; + k''pP)Q^' + d^pP'', (22) 
we combine the pure V7^-loop graphs, 

[(la) + (le)]o = +d^{ya{l^+mPLP''P'''] 

+pRya{H+mp''^m2r' +m+mPLP''^iPQ2r'- (23) 

The last two terms in the above are summed with the remaining two graphs to yield 
[(la) + (l.) + (lc) + (lJ)]o = +d^{Ya{Hi>)ypPLP''P'''] 

+PRya[{Hm+m)-ml)\p^^{PQ2)-' 

+ [(^ + m,-)(^+/^) -ml]ypPLP^P{PQ2)-\ (24) 

Since the above expression is sandwiched between the spinors of the initial and final states, it is 
tempting to replace m, in the last two terms by ^. But this is not legitimate as emphasized in the 
last section. Instead, m, should be replaced by ±^/2) on the rightmost (leftmost), in terms 
of the exact EoMs Q: 

u{p^)[{la) + (le) + (Ic) + (lJ)]o"(p+) 
= u{p^){d^{ya{li+i})ypPLP-'P''']+2ypPLP^PQ^')u{p^) 

+u{p-){lcd)qu{p+), (25) 

where the last term linear in q does not contribute to the charge but may contribute to the 
magnetic moment, 

{led), = lp^7^(^+/)^p^^(P22)-^-^^(^+/^)rp^LP^^(/'e2)-^ (26) 

In summary, leaving aside the {lcd)q term, we have 

/ 

£(U)o = +d^i{-m^^[{l^+^){mjPR + mlPL)-miml]{PQ2y^ 

x=a 

+ya{H^)ypPLP''''P''} +2ypPLP^Q2'- (27) 

The total derivative can be dropped in regularized loop integrals, so that Fig. 1 contributes to 
the Fi (0) term in Eq. ([T]) the following: 

+ eglY.\^ai\^ t YpPlP^Qi' = "regl I ypPLP^Q2\ (28) 

Qt Jk Jk 

where unitarity of V is used to finish the sum as the integrand is independent of m^. 
Now we manipulate iHfivi^)- First of all, {2b) and {2d) form a total derivative: 

{2b) + {2d) = -{cl-slr)^^,{kyQ2'). (29) 
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Using the shortcuts in Eqs. (I18I19I I. we have 

(2a) = 2cl,[5^{k^g^y-k^ky){QiQ2)-'-gMd-l)Q^'], (30) 
(2c) = 2cl;[^w{k^gnv-k^ky){QiQ2)-'+2{d-l)k^ky{Qi)~^]. (31) 

The last terms in (2a) and (2c) already form a total derivative. In the first terms, we decompose 
k^iQiQi)^^ = ^2 ^ +fnw(QiQ2)^^' then sum judiciously with (2/) to arrive at the result 

(2a) + (2c) + (2/) = 2ci{-{d-l)d^ikyQY')+g^y[Q^'+ml{QiQ2)-'] 

+k^ky [Q2^-{QiQ2)-'] -2k^kyQ2^^ 
= 2c^w{d^[kyQ2'-{d-l){kyQY')]+mlQ2'P^yy (32) 

Thus, using c^ +5^ = 1, the sum of all graphs is 

/ 

£(2x) = d^{kyQ2'-2cl{d-l){kyQ-')}+2mlQ2'P^y. (33) 

x=a 

Dropping the regularized total derivative and using Eqs. <^\5}, its contribution to the Fi (0) term 
in Eq. ([B is as follows, 

-egl I Q^'P^vfPL. (34) 
which cancels Eq. (|28|) . The vanishing charge is thus established at one loop in gauge. 



3.2 Magnetic moment 

Moving to the magnetic moment, we follow the computational procedure proposed in sec |2] 
Now only the graphs in Fig. 1 contribute. Since {lb) is quadratic in q when expanding in q, it 
drops out. The next simplest is (1/). Taking a derivative with respect to <7^, setting ^ = and 
making it manifestly antisymmetric in ix and v (denoted by the pair of square brackets below), 
we have 

[^v(l/)o] = +^^(«v+<v)('«'^i^ + '4^L)-2m,m2[r^,rv]), (35) 



where 



K% = i^[rAi,rv] + [rM,rv]A (36) 



K^y = m^.Yv] + [yii.yv]ii. (37) 



Anticipating that [^v(l/)o] is to be sandwiched between the initial and final spinors and noting 
that the 1^ in Kj^^ will yield a ^ upon loop integration, we can apply the limiting EoMs ]/)u — mtu 
after the moment has been isolated. The above is thus reduced to 

[^v'(l/)o] - +^^^^{K],y{m^ + ml)+2m,%,Yy]{m}-ml)), (38) 
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where, from now on, ^ means equality when sandwiched between the spinors or under the 
loop integration or both. All factors of Pm are removed in a similar fashion, confirming that 
the electric dipole moment does not arise at the one loop. 

Figs. (Ic) and {Id) should be treated together for symmetry reasons. There are two sources 
of terms, one from those explicitly linear in q and the other from the remaining terms (26) when 
computing the charge. Putting them together and taking the derivative, we have 

+\[f{^+^)f -f{li+^)f]PLgpvP^a{PQ2)~\ (39) 

where 

Pv;afi = -\dvPafi. (40) 
Antisymmetrization and applying EoMs yield, after some algebra, 

\d^{\c+\d)Q\ ^ niiKlJ — L — + 

-(ifi,+2m,[7^,rv])^^^, (41) 

where 

4v = h[lyv\+kv%M- (42) 

In deriving the above result, we used identities such as 

ynHyv-jvHJii = -\{H[yii.yv] + [y^.yv¥). (43) 
Myv-yvU = +\myv] + []^.yv]^). m 

Now we manipulate (le). Taking the derivative, plugging in the propagator PP'^ and doing 
antisymmetrization, one obtains 

[d^{le)o] = -i4v+4v)PL:^^+E,vPL^^^^. (45) 
where, using = m? and the identity 

to,7v]^ = k^[r^,rv]~2{[r^,^ky + [^,rv]k^), (46) 

the second term is recast as follows: 

E^y = i^{{HMyn.yv] + [yn.yvW+P))}f^ 

+2 { k^ {i> % yv] + %yvm+k,{^[y^M + [yM)]- (47) 
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Application of the limiting EoMs gives finally 

[dUle)o] - -(4, + 2m,-[7^,7v])^^ 

+ (^4, + {ml - mj)Kl, - 2mik^[r^,rv] +4m,-4v) 3^2^^^^ • (48) 



The graph Fig. 1(a) involves the triple gauge coupling and double gauge boson propagators, 
making it the most complicated to evaluate. We outline how this is accomplished. Taking the 
derivative and doing antisymmetrization we have 



where 



[<9v"(la)o] = ^P-^f^^+MPL^pa.^v. (49) 



ApcTTxv = A,hA^^^-^%^-h%^-\^^^^ (50) 



and 



^\v,pa — gvagfip — gvpgfia, (51) 
^Jiv,pa ~ k^{§vpka ~ gvakp) — kv[giipk(y — g^fykp). (52) 

The contraction with is standardized using Eq. (|43T ) into {K'^y +Kjj^y)Pi, while the contrac- 
tion with yields, by making use of Eq. (|44|) . 



{kM^, Tv] + Yv]^) + kvi^M + [rnm)PL, (53) 
which reduces to miK^y using EoMs. The final form is 

hk" 3 



[<9v'(la)o] - -{Kly + lmi[y^,Yv])[j^^^ + 



, ( hk^ h \ 
+lmiKly — + % — ■ (54) 

To summarize our calculation thus far, the terms relevant to the neutrino magnetic moment 
are given in Eqs. (|54I41I48I38I ). The next task is to demonstrate the E,yv cancellation among 
those terms. We first decompose 5^k^ = Qi — ^wQi to remove k^ from numerators in Eqs. 
(I54I48I) . The KJ^^ terms sum to 



I 2 f ^ 1 1 \ 



For any of the three terms in the above, the k^ and ky factors in Kjj^y may be simultaneously 
replaced by (k + p)^ and {k + p)y, because the resulted additional terms, upon the loop integra- 
tion, will be proportional to 

p^iy>,rv]+pv[YnJ], 
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which vanishes when sandwiched between the initial and final spinors. We make this replace- 
ment for the last term in the sum. Using again d^Qj^ = —Ik^Qj^ and d^P^ = —2{k + 
p)^P~^, the sum becomes 



The apparently <^vk -dependent terms in the sum ^ [(9y (lx)o], including the one in Eq. (1551) but 

x=a 

dropping total derivatives, are collected below: 

^ 5;: 



+ 



m,- 



m 



w 



[Kj^^m + [Tm. Tv] {m} - ml + ^wml)) . (56) 



The integral of the above second term is simplified using Eq. (|63l) and EoMs, while the first one 
is split by {QiQi)^^ = m^{Q^^ — ^2 ^^^^ dependence disappears completely 

from the sum: 



/ 



(57) 



PQi P^J' 

Adding the above with the terms that are explicitly i^vK -independent, we obtain the final sum of 
terms contributing to the neutrino magnetic moment: 

iwiM ^ "•.-te.rvlj^(^-^)-Wv+2".,te,rv])5/ 



+ ([m^^^(m^ -m7 



1 



(58) 



2]i^;v-6m,[7^,rv])g^- 

From Eqs. (II 1415 81) and the loop integrals defined in the appendix, we obtain for the neutrino 
V, the magnetic form factor at the vanishing momentum transfer. 



ai I 



a 



(59) 



where /i, 7i 2, ^1,2 are functions of the mass ratios Xa = in^/rn^ and yi = mj/in^. This result 
is indeed manifestly gauge independent in the class of gauges. 



4 Evaluation in unitary gauge 

Working in unitary gauge means that the limit i^iy — °° is taken before the loop integrals are 
evaluated. Since appears exclusively in the propagators of the gauge bosons, would-be 
Goldstone bosons and the ghosts c^, only the gauge boson propagator survives the limit, 

Pnv{k) ^ P^,y{k) = {g^v-m^%ky)Q^\ (60) 
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and thus only the pure-V7^ graphs (a, e) in Fig. 1 and (a, c) in Fig. 2 remain. We have 
presented our calculation in gauge in a way that can be easily adapted for unitary gauge. 

For the charge contribution from Fig. 1 we take the limit i^iy — > oo in the integrand (1231) 
where only the total derivative term survives: 

[(la) + (le)]o ^ d^{r^{}^+^)YpPLP-'pP''}, (61) 

whose integral vanishes in dimensional regularization. The relevant terms from the photon-Z 
mixing energy are obtained from Eqs. (13013 It . or more readily from Eq.(|33]). 

[(2a) + (2c)] ^ d^{-2ci{d-l){kyQY')}, (62) 

whose integral again vanishes. Thus the vanishing of charge at one loop occurs in unitary gauge 
in a stronger manner: each of the contributions from the proper vertices and the mixing energy 
vanishes separately. 

The magnetic form factor ^2(0) can also be obtained from intermediate steps in subsec |3^ 
We can sum Eqs. (I54I48I) and take the limit or cope directly with the total of all graphs 

since we know only Figs. 1 (a, e) survive the limit. The latter point can also be seen from explicit 
results in Eqs. (l41l38l) . Dropping the total derivatives and sending -t- 0°, the potentially ^ly- 
dependent part of the total in gauge, Eq. (l56l) . goes exactly to Eq. (l57l) without additional 
manipulations. The result in Eq. (|59| ) is thus recovered in unitary gauge. 

5 Summary 

The electromagnetic properties of neutrinos are an interesting topic that is potentially relevant to 
various astrophysical phenomena and laboratory measurements. Although we know from prin- 
ciples that the charge and dipole moments of a Dirac neutrino are physical quantities and cannot 
depend on computational methods or the choice of gauge in a consistent theory, this has never 
been explicitly examined before in a satisfactory manner even at one loop. We have studied this 
issue in the minimally extended standard model that incorporates neutrinos masses and mixing. 
We demonstrated at one loop in both R^ and unitary gauges that the magnetic moment and 
vanishing charge are indeed gauge-independent quantities. This statement is exact in the sense 
that it is true for any values of various masses and the lepton mixing matrix as long as the latter 
is unitary. We have accomplished this by manipulating directly the integrands of loop integrals 
and employing simple algebraic identities like (I43I44I46I) and integral relations like (|63T ). We 
believe this approach is advantageous over the one that handles the results of loop integration, 
and may be useful in other contexts. Finally, we mention that various approximations to our 
exact one-loop result for the magnetic moment in Eq. (1591 ) are possible. For instance, when all 
neutrinos and charged leptons are much lighter than the weak gauge bosons as is the case in SM, 
we have from the explicit results in the appendix that F2(0) ~ —{SGFtnj) / {4-n^\/2), where the 
mixing matrix drops out from the leading term, so that the interaction potential of the neutrino 
V,- of mass m,- and spin S with an external magnetic moment is, V ^ {SeGfirii) / {4n^\/2)S ■ B, 
recovering the well-known result in the literature. 
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Appendix: some useful integrals 

We list some loop integrals relevant to our evaluation of the magnetic moment. The follow- 
ing relation is used in sec |3]for reduction of terms: 

^ 2 [ kg ^ 1 m\-ml-p^ \ 

where D\ = {k-\-pY — m\, D2 = — m\. Using notations in Eq. (|20|) with — mf, the basic 
integral is 

where Xg = m^/m^ and yi = mj/myy. For simplicity, we also define the integrals 

1 1 \ i 



PQ, - i4n/'^'"^''^^ ^^^^ 

1 z 1 , , 

■/lUa,y/), (66) 



kP^Qi {^Ttfml 

^^(^t'^^^^"'^'^' '''' 

2k^ _ i 



\ PQ\ (4;r)2 -2 



m 



K2{xg,yi). (69) 



w 



The parametric integral for I{s,t) is 

I{x,y) = J dt \n[xt + {\-t)-yt{\-t)-iO^]. 
The other functions are related to it by 



(70) 



Ii{x,y) = lnx-I{x,y), (71) 

Jiix,y) = j^I{^,y), (72) 

J2{x,y) = Ji{l/x,y/x), (73) 

Kiix.y) = y-'[hix,y) + {l+y-x)Ji{x,y)l (74) 

K2{x,y) = y~'[I{x,y) + {l+y-x)J2{x,y)]. (75) 
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Note that the singularity at 3; = is spurious since the original integrals are smooth there. 

The analytic result for / is known for all parameter regions, but we only record it for the 
case relevant to SM, i.e., for < y < x <^ 1, 

I{x,y) = -2-^{l-x-y)lnx+^lnR, (76) 
2y 2y 

where 

X = {l+x'W~2x-2y-2xyy/\R= . (77) 

1 +x — y + A, 

The other two functions are 

Ji{x,y) = -^^^\nR + ^\nx, (78) 
2yA 2y 

JiM = -Ji{x,y)-j\nR. (79) 
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